Inelastic neutron scattering studies of Crystal Field Levels in PrOs 4 Asi 2 
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We use neutron scattering to study the Pr 3+ crystalline electric field (CEF) excitations in the 
filled skutterudite PrOs4Asi2. By comparing the observed levels and their strengths under neutron 
excitation with the theoretical spectrum and neutron excitation intensities, we identify the Pr 3+ 
CEF levels, and show that the ground state is a magnetic T 4 triplet, and the excited states hi, 
I4 and T23 are at 0.4, 13 and 23 meV, respectively. A comparison of the observed CEF levels in 
PrOs4Asi2 with the heavy fermion superconductor PrOs4Sbi2 reveals the microscopic origin of the 
differences in the ground states of these two filled skutterudites. 

PACS numbers: 75.47.-m, 71.70.Ch 



I. INTRODUCTION 

The Pr-based filled skutterudites (FS) have the for- 
mula PrT 4 Xi2, where T is one of the transition met- 
als Fe, Ru, or Os, and X is a pnictogen (P, As, or 
Sb)^2i£. The notably mounting interests and efforts in 
the study of the FS compounds are motivated by the 
remarkable diversity of their electronic and magnetic 
ground states, including multipole ordering 4 ^, small 
gap insulators^' 7 , conventional superconductivity^, 
unconventional superconductivit y 10 ' 11 and magnetic 
orderin g 12 : 13 i 14 : 15 i 16 . Despite the large differences in 
their physical properties, these compounds are gov- 
erned by only a few parameters, including the inter- 
action between the conduction and the 4/ shell elec- 
trons (the c-f coupling) and the effect of the crys- 
talline electric field (CEF) potential on the Pr a+ if 
electronsi^^^ii 1 .^^^^^ 7 . For example, trans- 
port and bulk magnetic measurements on the heavy 
Fermion superconductor PrOs4Sbi2 suggested either a Ti 
singlet ground state or a nonmagnetic doublet ground 
stat o 10 ' 11 . Inelastic neutron scattering experiments on 
PrOs 4 Sbi2 showed that the Pr 3+ CEF levels include a 

(2) 

Ti singlet ground state and a low-lying Y\ magnetic 
triplet excited state at 0.6 me V 18 ' 19 ' 20 . This rules out 
the quadrupolar Kondo effect, which arises only from a 
nonmagnetic doublet ground stat&i 7 -, as the microscopic 
origin for the observed heavy-fermion superconductivity. 

The FS compounds belong to the space group Im3 3 . 
The rare earth atoms are located at the corners and 
body-center of the cubic lattice, each of which is sur- 
rounded by a simple cube of 8 transition metal atoms at 
the 8c sites [Fig. 1(a)] and by an icosahedron of 12 pnic- 
togen atoms at the 24g Wyckoff sites [Fig. 1(c)]. Owing 
to their unique structure, a subtle modification on com- 
position can result in a different CEF scheme and thus 



a completely different ground state. However, a general 
understanding is desirable as to how the compositions in- 
fluece the CEF levels. In PrOs4Asi2, in which the pnicto- 
gen Sb in PrOs 4 Sbi2 is replaced by As, the material dis- 
plays quite different correlated electron propertie s 14 ' 15 . 
The temperature dependence of the electrical resistiv- 
ity reveals Kondo lattice behavior, which is not ob- 
served in PrOs 4 Sbi2 21 . Specific heat measurements in- 
dicate an enhanced electronic specific heat coefficient of 
7 w 1 J/mol K 2 for T < 1.6 K and < H < 1.25 Til 
The compound exhibits several ordered phases at tem- 
peratures below 2.3 K and fields below about 3 T— . The 
ground state has been determined to be antiferromag- 
netic (AF) by neutron scattering experiments 15 . A de- 
termination of the Pr 3+ CEF level scheme in PrOs 4 Asi2 
and its microscopic origin is crucial for understanding 
why its ground state is different from that in PrOs 4 Sbi2 . 
The outcome will lead to a more general understanding 
of how the structures and compositions in Pr-based FSs 
can influence their CEF levels and ground states. 

II. EXPERIMENTAL 

PrOs 4 Asi2 single crystals were grown using the molten 
metal flux method described in ReL 1 ^ and crushed 
into fine powder. Our neutron scattering experiments 
were carried out on the cold neutron triple-axis spec- 
trometer SPINS at the NIST Center for Neutron Re- 
search (NCNR) and on the HET chopper spectrome- 
ter at ISIS (Rutherford Appleton Laboratory), as de- 
scribed previously^ 3 .. We reference positions in recip- 
rocal space at wave vector Q = (q x ,QyiQz) in A -1 us- 
ing (H,K,L) reciprocal lattice units (r.l.u.) notation, 
where (H,K,L) = (q x a/2ir 1 q y a/2ir, q z a/2ir) for the cu- 
bic PrOs 4 Asi2 unit cell (a = 8.5319 A)il We used a 
3 He- 4 He dilution refrigerator for the field-dependent ex- 
periments. The nature of observed CEF excitations were 
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FIG. 1: (Color online) (a) The cube of 8 Os ions surround- 
ing the central Pr 3+ ion in PrOs4Asi2- These give an Oh- 
symmetric CEF. (b) The corresponding spectrum of Oh- 
symmetry Pr 3+ CEF levels, (black =singlet, blue=doublet, 
red=triplet). The relative coupling x that gives singlet-triplet 
degeneracy is shown by a dashed vertical, (c) The 12 nearest- 
neighbor As ions surrounding the central Pr 3+ in PrOs4 Asi2 , 
giving a reduced symmetry (Th) CEF. (d) The corresponding 
As-only TVsymmetry Pr 3+ CEF spectrum in PrOs4Asi2. 



confirmed in a large temperature (0.08 K-200 K) and 
magnetic field (0 T-ll T) range. 

III. RESULTS AND DISCUSSION 

Figure 2 summarizes the neutron scattering intensity 
from PrOs4Asi2 on HET at temperatures between 1.5 K 
and 200 K. Since the CEF magnetic scattering decreases 
with increasing Q whereas the intensity of phonons in- 
creases with Q, a comparison of the neutron intensities 
in the low- and high-angle detectors can distinguish be- 
tween magnetic and phonon scattering. Figure 2a shows 
the scattering function at T — 1.5 K and T = 200 K with 
an incident neutron beam energy of Ei = 32 meV. Com- 
parison of the low- and high-angle data reveals two clear 
CEF excitations at 13 meV and 23 meV, with phonons 
at ~ 20 meV. Measurements with Ei = 12 and 50 meV 
showed no evidence of additional CEF excitations at en- 
ergy transfers between 2 and 8 meV or above 25 meV 
[Figs. 2(b) and 2(c)]. 

To search for CEF excitations at energies below 
2 meV, we carried out high resolution measurements us- 
ing SPINS. At T = 0.32 K, energy scans at Q = (1.2, 0, 0) 
showed a clear peak at 0.4 meV; this mode decreases and 
becomes broader on warming to 2.5 K and 6 K [Fig. 3(a)] . 
Figure 3(b) shows that the energy of the ~0.4 meV mode 
has weak Q-dependence and decreases in intensity with 
increasing Q, thus confirming its magnetic nature. Fig- 
ure 3(d) reveals that the elastic intensity also decreases 
on warming from 0.08 K to 4 K. This reduction of inten- 
sity in the elastic channel with increasing temperature is 
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FIG. 2: (Color online) (a) Neutron inelastic scattering at 
T — 1.5 K and T = 200 K with E z = 32 meV, integrated 
over scattering angles from 9° to 30° (low-angle detectors) 
and from 105° to 140° (high-angle detectors), (b) The same 
at Ei = 12 meV; (c) Ei = 50 meV. The scattering function 
S(Q,u) was normalized by comparison to a vanadium stan- 
dard. 



also observed in the HET data, evidencing the ground 
state is a magnetic multiplet. 

Figures 4(a)-(d) show the temperature dependence of 
the low-angle scattering for Ei = 32 meV. The CEF peak 
intensities do not change significantly with temperature 
between 1.5 K and 5 K. At 50 K the intensity in the elas- 
tic channel has undergone a substantial decrease, and the 
13 meV peak has shifted to 10 meV. On further increas- 
ing the temperature to 100 K and 200 K the intensities 
at meV, 13 meV and 23 meV continue to decrease, 
whereas the scattering at 10 meV increases. 

The theoretical description of the Pr 3+ CEF levels in 
PrOs4Asi2 is complicated by the presence of important 
contributions from two sets of neighboring ions, Os and 
As. The Pr 3+ ion in Pr-based FS has a 4/ 2 configuration, 
which in Russell-Saunders coupling has a ninefold degen- 
erate 3 H4 ground state. This degeneracy is lifted by the 
CEF interaction, which we assume to be dominated by 
the 12 nearest neighbor pnictogens (As) and the 8 next 
nearest neighbor Os ions; the distances to these ions are 
dpr-As = 3.23 A and dp r -Os = 3.69 A respectively. 

A. Single — charge model with separate ions 

Os ions form a simple cube around the Pr ions and 
they alone give an Oh symmetric CEF. The arrangement 
of the 12 pnictogens (As) around Pr 3+ forms 3 orthogo- 
nally intersected planes where the As- As bonds are shown 
as solid lines with length L and W is the length of the 
dashed lines in Fig. 1(c) {b = W/L = 0.4267 ^ 1). When 
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the 4 pnictogen atoms in each of the 3 orthogonally in- 
tersected planes form a square, i.e., b = W/L — 1, the 
fourfold rotational symmetry is recovered and the point 
group symmetry becomes Oh with the simple cubic CEF 
potentials [Fig. 1(c) This Oh case is treated by Lea, 
Leask, and Wolf (LLW) 26 (see their Fig. 9); we have red- 
erived their excitation spectrum as shown in Fig. 1(b). 

Both O h and T h CEF interactions split the Pr 3+ 3 H 4 
ground state into a singlet, a doublet and two triplets^ 5 .. 
In the 6 =/= 1 T/j-symmetry case, these multiplets are 
referred to as T\ (a singlet), T 2 3 (a nonmagnetic dou- 
blet), and I4 and I4 (magnetic triplets). These 
two triplets are linear combinations of the O^-symmetry 
triplets, mixed by the new Th CEF interaction^ 5 -; this 
mixing modifies the excitation spectrum and leads to 6- 
dependent neutron transition intensities. 

The X/j-symmetry CEF excitation spectrum has not 
been considered in detail in the literature, and the corre- 
sponding neutron transition intensities between Th CEF 
levels have not been considered at all. To aid in the inter- 
pretation of our neutron scattering data we carried out 
these CEF calculations using a point charge model. We 
assumed an expansion of the perturbing CEF potential 
in spherical harmonics, 

e 

V(Sl) =J2^Y1 M im Y em (fl), (1) 

1=4,6 m=-l 

where the interaction strengths 54 and g 6 are treated 
as free parameters. The spherical harmonic moments 
{Mi m } are determined by the positions of the 12 As 
ions, which we assigned the (scaled) coordinates x = 
(±1,±6,0), (0,±1,±6), (±6,0, ±1). The nonzero inde- 
pendent moments for I = 4, 6 are M40 = 21(1 — 3b 2 + 
6 4 )/20F(l + 6 2 ) 2 , M 60 = 3\/l3(2- 176 2 + 26 4 )/80F(l + 
6 2 ) 2 and A/ 66 = -15%/30W(l-6 2 )/16^F(l + 6 2 ) 3 . The 
nonzero Mqq for b 7^ 1 (Th symmetry) confirms the pres- 
ence of the Bf terms of Takegahara et al. [Eq.(7) of 
Refi^S,] , in addition to the usual B 4 and B\ O^-symmetry 
terms. (Note that the TVallowed moment Mqq vanishes 
at the Oft,-symmetry point 6=1.) We also confirmed that 
the other nonzero moments satisfy the ratios quoted in 
Eq.(7) of Ref. 25 . Unlike Takegahara et al£i, we do not 
introduce a new parameter y for the X/j-symmetry terms, 
because they are completely determined by g§ and the 
lattice parameter 6 in the point charge model. This was 
previously noted by Goremychkin et al*&. 

Diagonalization of this Th CEF interaction within the 
p r +3 3 jj 4 nonet gives our results for the spectrum of CEF 
levels and their associated eigenvectors. These eigen- 
vectors depend only on the ratio g^/gi and the lattice 
parameter 6; the energies in addition have an arbitrary 
overall scale. Our results for the spectrum for 6=1 
(Oh symmetry) and 6 = 0.4267 (PrOs4Asi2 geometry) 
are shown in Figs. 1(b) and 1(d), using LLW normal- 
ization conventions^ 6 -. (These conventions set our two 
Hamiltonian parameters in Eq.(l) to 54 = (9687r/21)a; 
and 5 6 = (-5808tt/221)(1 - \x\).) Note that the 6 = 1 



and 6 = 0.4267 level schemes differ qualitatively, which 
demonstrates the importance of the Th terms in this 
problem. 

We find that the Oh singlet and doublet energy eigen- 
vectors are unmodified by the Th interaction, consistent 
with Takegahara et al£^. The singlet eigenvector (in a 
Jf 4 basis) is |*i) = V^7l2|0) + V5/24(\4) + | - 4)) 
and the two doublet states are | v E , 23a) = — \/5/12|0) + 
77/24(14) + |- 4)) and |* 236 ) = VT^G 2 ) + I - 2)), con- 
sistent with earlier (numerical) result s 25 ' 26 . The singlet 
and doublet energy eigenvalues in our conventions are 
modified by the Th interaction. In terms of the LLW vari- 
able as^ 6 - and our parameter 6 they are Ei = —(16/13(1 + 
6 2 ) 2 )(91x(l - 36 2 + 6 4 ) - 20(1 - |x|)(2 - 176 2 + 26 4 )) and 
£23 - -(16/13(l + 6 2 ) 2 )(13x(l-36 2 + 6 4 ) + 16(l-|a;|)(2- 
176 2 + 26 4 )). The corresponding analytic results for the 
two Th triplet states for general 6 are quite complicated, 
so we only present numerical results for these states. 

The neutron transition intensities are defined by X/ — 
K/I«^°'l*>| 2 > as introduced by Bir geneav*2£. (There is 
an implicit sum over initial and final magnetic quantum 
numbers.) Our X^-symmetry results for these quantities 
are shown in Fig. 5(a). The values in the limits x = ±1 
(no £ = 6 term, hence Oh symmetry) implicitly check Bir- 
geneau's numerical Oh results; see the off-diagonal entries 
in his Table 1(e). These Oh limits are indicated on the 
vertical axis of Fig. 5(a). 

Next we compare the observed CEF levels and their 
neutron excitation intensities to the well-known LLW 
CEF results for Oh symmetry [Fig. 1(b)] and our calcu- 
lated CEF predictions for PrOs 4 Asi2 under Th symmetry 
[Fig. 1(d) and Fig. 5(a)]. Both O h and T h CEF spectra 
have x values that can accommodate a magnetic triplet 
ground state and a nearly degenerate singlet first excited 
state [vertical lines in Figs. 1(b) and 1(d)]. However, it 

is evident that the Oh scheme cannot explain the data 

(2) 

because the observed 0.4 meV transition T\ J — ► Tj is in- 
correctly predicted to have zero intensity due to the Oh 
symmetry. While in the Th scheme the relative neutron 
excitation strengths of the higher levels (at 13 and 23 
meV) predicted in Fig. 1(d) seem to be in good agree- 
ment with observation at low temperatures, the As CEF 
alone predicts an incorrect spectrum of levels [Fig. 1(d)], 
with the triplet I4 being the highest excitation. The 
calculated neutron transition intensity shown in Fig. 5(a) 
can not explain the observed intensity at higher temper- 
atures. As temperature increases, the excited states get 
populated and the excitations start to decrease in inten- 
sities. Meanwhile the new excited-state transitions start 
to increase. If instead of T23, is the highest level, 
the intensity at 12.6 meV would not increase but that 
at 22.6 meV would, because the Ti to T 2 3 transition is 
not allowed even in Th symmetry. Goremychkin et al^ 
showed that the highest level in the similar Sb material 
is the r 2 3 doublet. 

B. Combined Os As CEF model 
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FIG. 3: (Color online) a) Low energy spectrum of CEF ex- 
citations observed at T = 0.32, 2.5 and 6.0 K using the 
SPINS spectrometer at NCNR. (b) The wave vector depen- 
dence of the excitations at Q = (0.8, 0, 0), Q = (1.2, 0, 0), and 
Q = (1.7,0,0). (c) The expected and observed temperature 
dependence of the intensity of the 0.4 meV mode, (d) The 
temperature difference spectrum between 0.08 K and 4 K, 
showing clear reduction in magnetic elastic scattering. 
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FIG. 4: (Color online) The temperature dependence of the 
excitations observed on HET with Ei — 32 meV at (a) T = 
1.5 K, 5 K; (b) 50 K; (c) 100 K, and (d) 200 K. The lines are 
theoretical results for neutron excitation intensities, from the 
combined Os-As CEF model, with an arbitrary overall scale 
factor. 



The twin constraints of having the 1^3 level at the top 

(2) 

of the spectrum and having a large T 4 «-> Fi neutron 
excitation strength requires both Os and As terms in the 
CEF interaction. We therefore introduce a combined Os- 
As Hamiltonian, 



H = H(Os)+H(As). 



(2) 



Although this model nominally has four parameters 
(<7. 
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FIG. 5: (a) The theoretical neutron transition intensity for As 
CEF alone with b=0.4267. (b) Neutron excitation intensities 
predicted by the combined Os-As CEF model. 



g° s and gf s cannot be distinguished because they are 
summed into a single coefficient of the O/j-symmetry 
I = 4 interaction. For this reason we introduce com- 
bined O/i-symmetry Os-As coefficients 34 



n Os 
94 



As 
9i 



28.g 4 - out/6, 
- 64g6, where 



and <?6 = g@ s + g@ S : which we normalize according to 
LLW conventions. As the I = 6 X/j-symmetry terms 
from H(As) in the CEF are proportional to gif s alone, 
the strengths g° s and gif s can be distinguished. We 
parametrize these two £ — 6 interactions using the to- 
tal O/j-symmetry g 6 and a Th/Oh relative strength 
which is the ratio of the coefficients of 162 to Yqq in 
the CEF potential. The energy levels of this Hamilto- 
nian are E^) = -§9i ~ &96 - f, Ei = 
#4(1) = -6ff4 - 8g 6 + / and E23 = 4g 4 
/ = ((20g 4 + 12g 6 ) 2 + 960 r^f) 1 / 2 . For r 6 = these 
reduce to the familiar LLW Oh spectrum. In the pure 
As model, r§ is determined by CEF theory if we assume 
point As ions, and is given by (11\/105/4)6 2 (1 — b 2 )/{\ + 
fe 2 )(l-(17/2)& 2 +6 4 ). For PrOs 4 Asi2 we have b = 0.4267, 
which gives a rather large r§ — —6.901. This drives 
strong level repulsion between the two triplets, which 
explains why the pure As spectrum of Fig. Id differs so 
greatly from the Oh (pure Os) symmetry spectrum of 
Fig.lb. 

Our experimentally observed CEF levels are close to 
but not exactly consistent with the predictions above of 
the mixed Os-As model, since the gap ratio (£4(1) — 
#4(2))/ (#23 — #4(2)) ~ 0.57 is slightly below the theo- 
retical lower bound of 7/12. The parameters we esti- 
mate from the measured gaps are 54 « 0.24 meV and 
ge « 0.20 meV. The value of tq is not determined by the 
measured energies due to the inconsistency mentioned 
above, although < 0.5 appears plausible. A more sen- 
sitive determination of r§ is possible through the mea- 
surement of the inelastic neutron excitation intensities 
we discuss below. 

The neutron excitation intensities in this combined Os- 
As Hamiltonian depend only on a single parameter 9, 
which is the mixing angle of the triplet energy eigenvec- 
tors when expanded in an O^-symmetry |3), |3') basis, 



|4(1)) = +sin(0)|3) + cos(0)|3') 
1-4(2)} = + cos(0)|3) -sin(0)|3') 



(3) 



Os 



n Os 
96 ■ 



As 
54 



and <?g ), only three are independent; This mixing angle is related to the Hamiltonian param- 
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eters by tan(20) = 2\Zl5r 6 /(5(g 4 /.9 6 ) + 3). The sin- 
glet and doublet Oh energy eigenvectors are unchanged. 
The nonzero neutron excitation intensities in terms of 
s = sin(0) and c = cos(0) are r£ 2) T x = (20/3)s 2 , 

rf <-» r« = 7/2 + 8cV, r£ 2) «-» r 23 = 4 + (ie/3) s 2 , 
ri «-» rf } = (20/3)c 2 , r 4 1} «-» r 23 = 28/3 - (16/3)* 2 , 

rf «-> rf = (25/2) (1 - (4/5)s 2 ) 2 , and v[ 1] «-> r£ x) = 
(1/2)(1 + 4s 2 ) 2 . The calculated neutron scattering inten- 
sity of different transitions as a function of is shown 
in Fig. 5(b). We recover the O^-symmetry results of 
Birgeneau (Table 1(e) of Refill) for s = 0, c = 1. 

We carried out a least-squares fit of our neutron exci- 
tation data at 1.5 K, 50 K, 100 K and 200 K (Fig. 4) 
to the theoretical intensities given above, which gives an 
estimate of the triplet mixing angle in PrOs4Asi2, 

0w22.5°. (4) 

When combined with the values of 54 and g§ from the 
spectrum, this corresponds to rg « 1.2. In this fit 
the relatively isolated rjj 2 — > T^p peak at 23 meV was 
used to infer the background, which was taken to be con- 
stant plus linear. The assumed lineshapes were Loren- 
zians with a common lincwidth, fixed by the 23 meV 
peak. The calculated intensities of the individual transi- 
tions (dotted lines) and the total intensity (solid lines) for 
each temperature are shown in Figs. 4(a)-(d). We note 
that the intensity reduction at 0.4 meV on warming from 
0.32 K to 2.5 K is larger than that expected from the CEF 
model [Fig. 3(c)], thus suggesting Pr-Pr interactions be- 
low Tjy (=2.3 K) are important. On the other hand, the 
large difference between the calculated and expected in- 
tensity around 8 meV in the 200 K data is presumably 
due to thermally populated phonons (Fig. 2a). 

C. Field effect on the CEF gap 

(2) 

If the ground state of PrOs4Asi2 is indeed the P^ 
triplet, application of a magnetic field should Zccman 
split it, resulting in a field dependent energy gap. There 

should also be a reduction in the intensity of the zero- 

(2) (2) 

energy — > T\ magnetic scattering. Figure 6 shows 
that these expectations are indeed qualitatively satisfied. 
The first excited state at 0.4 meV shifts toward higher en- 
ergies as the applied field increases. The field-dependent 
transition energy is linear only at higher fields (between 
6 T and 11 T). The drop of intensity in the elastic chan- 
nel is almost constant for all applied fields, as shown 
in the inset of Fig. 6(b). Figure 6(c) shows that the 
wave vector dependence is also present with applied mag- 
netic field {H = 9 T). Normally the field-splitting of the 

ground state multiplet would be a very clear test of our 
(2) 

T\ triplet assignment for the ground state. However, 

PrOs4Asi2 is complicated by the near degeneracy of the 
(2) 

T\ and T\ levels, which mix strongly under an applied 
field. This results in a more complicated spectrum of 
low-lying states, with several low-field level crossings and 
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FIG. 6: (Color online) (a) The magnetic field dependence 
of the low-lying CEF excitations at Q = (1.2,0,0) and 
T — 0.08 K. The inset shows the field dependence of the first 
excited state, (b) The difference spectrum between T and 
11 T at Q = (1.2, 0, 0) and T = 0.08 K. The effect of an ap- 
plied field is to suppress intensity at hui — meV, and to split 
the spin-triplet ground state; the latter results in the field 
dependence of the 0.4 meV peak in (a), which may involve 
an intra-triplet transition. The elastic intensity suppression 
effect essentially disappears for fields above 4 T. (c) The field- 
induced CEF excitation at « 1.1 meV is weakly wave vector 
dependent, and shows essentially no temperature dependence 
between T = 0.08 K and 4 K. 
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FIG. 7: (Color online) The boundary between singlet and 
triplet ground states in skutterudites (Ei = E^) as a func- 
tion of re and Qi/g%, and the observed PrOs4Asi2 spectrum. 

neutron scattering intensities that are also modified by 

(2) 

their field-induced T\ -T\ mixing. 

Our determination of the CEF levels in PrOs4Asi2 
reveals the reasons for the wide range of behaviours in 
different FS. The spectrum of CEF levels is largely de- 
termined by the Oh symmetry field of the eight nearest 
neighbor ions, and for Os the near equality of the I = 4 
and £ — 6 strengths (74 and g§ implies nearly degener- 



6 



ate low-lying singlet (insulator) and triplet (AF) levels. 
The low temperature magnetic properties are determined 
by which of these phases happens to be the true ground 
state. In the CEF model, this is specified by the two 
parameters g^/ge and r§ (Fig. 7); in PrOs4Asi2, which 
has a triplet ground state, we estimate 54/^6 ~ 1-15 and 
rg w 1.2. The Th symmetry pnictogen CEF (propor- 
tional to re) acts to stabilize the triplet state, and can 
itself lead to a triplet ground state if r 6 is sufficiently 
large to cross the phase boundary shown in Fig. 7. 

In principle, one can extend our approach to calculate 
the ground states of other Pr-FS by determining its crys- 
tal structure and 34/36 ratio. The necessity of using the 
T/j-symmetry of As rather than the O/j-symmetry pure 
Os form to explain the observed excitations shows that 
the detailed pnictogen geometry is important in deter- 
mining the CEF levels. Indeed, the AF-ordered ground 

(2) 

state in PrOs4Asi2 can arise from a T\ triplet mag- 
netic ground state, while the superconducting PrOs4Sbi2 
has nonmagnetic Ti singlet ground state. The nearly 
degenerate first excited state Ti at 0.4 meV (~ 4 K), 
and its temperature and field dependence (Figs. 4 and 
6), may explain the presence of multiple transitions in 
the specific heat (in C(T)/T versus T) and its field 
dependencei^£i22. 

II. SUMMARY 

To understand the observed Pr 3+ CEF levels, one 



must incorporate the As ions' contribution to the CEF 
Hamiltonia n 24 ' 25 , in addition to the usual Os cubic field 
terms. A comparison of our CEF calculations using this 
more general Hamiltonian with our experimental results 

shows that the Pr 3+ CEF level scheme in PrOs4Asi2 con- 

(2) 

sists of a magnetic triplet ground state, a nearly 

degenerate T\ singlet excitation, and higher r^ 1 ' mag- 
netic triplet and T23 nonmagnetic doublet excited states. 
We find that contributions in the CEF Hamiltonian due 
to As are important in determining the neutron exci- 
tation intensities in PrOs4Asi2; our results differ qual- 
itatively from the predictions of the conventional CEF 
Hamiltonia n 26 ' 27 , and therefore provide a microscopic 
understanding for its AF ground state. 
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